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We present a generalized model of a diffusion-reaction system where the reaction occurs only 
on the boundary. This model reduces to that of Barato and Hinrichsen when the occupancy of 
the boundary site is restricted to zero or one. In the limit when there is no restriction on the 
occupancy of the boundary site, the model reduces to an age dependent Galton- Watson branching 
process and admits an analytic solution. The model displays a boundary-induced phase transition 
into an absorbing state with rational critical exponents and exhibits aging at criticality below a 
certain fractal dimension of the diffusion process. Surprisingly the behavior in the critical regime 
for intermediate occupancy restriction A'^ varies with N . In fact, by varying the lifetime of the active 
boundary particle or the diffusion coefficient in the bulk, the critical exponents can be continuously 
modified. 

PACS numbers: 64.60.Ht, 64.70.-p, 68.35.Rh 



Nonequilibrium phase transitions can differ signifi- 
cantly from their equilibrium cousins, and are a subject 
of continuing interest. Recently, Barato and Hinrich- 
sen [U [2] (BH) studied a reaction-diffusion model ap- 
parently exhibiting a new universality class of nonequi- 
librium phase transitions, boundary induced phase tran- 
sitions. The model they studied was a variant of a model 
introduced earlier by Deloubriere and van Wijland [3, 
who however did not discover the novel scaling behavior. 
In this letter we generalize the model and solve it ana- 
lytically in a certain limit. Our model exhibits a phase 
transition and at criticality we find the presence of ag- 
ing [H [S] , a nonequilibrium property observed in such 
diverse systems as spin-glasses gels [5] and turbu- 
lence [3 . Furthermore a critical fractional dimension for 
the diffusion process is naturally obtained in the con- 
text of our model. Interestingly we find that the critical 
exponents of the model vary continuously with the pa- 
rameters of the model. We also note that our toy model 
is interesting in the context of catalytic reactions and bi- 
ological situations when the reaction occurs on a specific 
membrane, such as a boundary of a cell. 

Our model is defined as follows: As in the BH model, 
a particle starting on the boundary follows a birth-death 
process. It produces an offspring (next to the bound- 
ary site) with probability p {A^A -I- O) and dies with 

probability 1 — p {A >0). The particle A continues to 

reproduce until it dies. The offspring O diffuse freely in 
the bulk since their birth and up to the time they reach 
the boundary, and then they start to reproduce again 
according to the same death-birth process. By varying 
the rate of offspring production, one expects to reach ex- 
tinction for p — >■ (absorbing state) and a growing or 
stable population for p — > 1, intuitively one expects a 
transition between those two states as p grows from 
to 1. In the BH model, there is a constraint that only 
one particle can exist on the boundary site at any give 
time. Any other particle trying to enter the boundary 
site dies immediately. We generalize this to allow up to 



N particles to coexist simultaneously on the boundary. 
As shown by BH, the role of the bulk diffusion is just 
to generate a probability distribution ip{t) for the arrival 
of the newly born particle to the boundary site, where 
it can, if the boundary site is not fully occupied, repro- 
duce some number of times, setting off new processes and 
then die. Introducing constraints on the occupancy of the 
bulk sites has no effect on the global dynamics, due to 
the indistinguishability of different particles [H]. By to- 
tally relaxing the constraint on the boundary occupancy, 
i.e., N — )■ 00, we are capable of completely solving the 
model since technically the unrestricted model is a ran- 
dom branching, or Galton- Watson, process. 

In the unrestricted case we can make another simplifi- 
cation, since the time to die is governed by a short-range 
exponential distribution, we can, without doing any harm 
to the model, consider all births as happening simultane- 
ously. In this view, the particle generates k children, gov- 
erned by a geometrical distribution, pk = (1 — p)p^ , and 
then immediately dies. The basic technique employed in 
the analysis is the use of the generating functional [9] 
and parallels the solution of the standard age-dependent 
branching process HUj. Let us denote by Z{t,t + s) the 
number of particles that been observed on the boundary 
during the time interval {t,t + s). We define a generating 
function for the random variable Z{t, t-\-s) as the average 
of s^(*'*+«) over the distribution of Z{t, t + s) 

OO 

G(,,,+,)(,s)^ ^ P(Z(i,t + s)).^(*'*+^), (1) 

Z(t,t+s)=0 

where P{Z{t,t + s)) denotes the probability to observe 
Z(t, t + s) particles on the boundary in the specified time 
interval. The generating function for the number of off- 
springs is Go{s) — (s'^) = (1 — p)/{l — ps). We set the 
initial conditions such that at t = a single particle was 
injected into the bulk and it has the probability distri- 
bution ipit) to return to the boundary at time t — t^. 
By conditioning on the outcome of the returning time 
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for the first particle we obtain [TI] the equation for 

G{t,t+s){s) 



G(t,t+s)(s) = / Go{G(^t-u.t-u+s){s))^{u)du 

"'0 

rt-\-s /'OO 

+s / G'o(G(o,t-«+s)(s))V'(")rfM + / ^{u)du. 

Jt Jt+s 

(2) 

In the limit of s — > the generating function goes to 
P{Z{t,t + s) = 0) := PE{t,s), i.e. the probability that 
not a single particle appears on the boundary at the men- 
tioned time interval, taking this limit in Eq. ([2| we obtain 



PEit,s) 



1 P /"oo (3) 

ip{u)du+ / ip{u)du, 



I - pPsit- u,s) 



t+s 



where we have now used the explicit from of Go(s). 
Eq©, a non-linear Volterra equation of the second kind, 
is our main equation for the unrestricted case of the 
model and it describes the time evolution of a two-time 
quantity. We haven't assumed anything as to the form 
of 'ip{t) and so Eq. ([S]) is quite general. The long-time 
behavior of Pe is, as we shall see, governed entirely by 
the long time behavior of tp{t). For normal diffusion in 
the bulk, this long-time behavior is tpt ^ tpoot^^'^^^K 
with /3 = 1/2, ijjoo = l/V47rD, where D is the diffusion 
constant (defined as usual). We also consider the more 
general case of < /3 < 1, and in such case the con- 
stant ipoo is used for the normalization of tpt- The case 
of < < 1 can be achieved if the diffusion in the bulk 
is anomalous [12]; e.g., subdiffusion for < /3 < 1/2 de- 
scribed by models such as the continuous time random 
walk (CTRW) [HdS]. Herein we study Pe exclusively; 
other quantities such as the mean number of particles 
can also be obtained from Eq. ^ [TT]. 

We will now explore the behavior of PE{t,s) in the 
asymptotic limit of large t for different values of the pa- 
rameters p and /3. Inspired by numerical solutions of Eq. 
([3]), we adopt the ansatz 



PE{t,s)^P^ -A{s)t- 



(4) 



In order to obtain a solution of Eq. ([s]) we substitute the 
asymptotic form in Eq. Q into both sides of Eq. ([S]), 
perform an expansion in inverse powers of t and com- 
pare the coefficients in front of the appropriate leading 
terms. We leave the exact technical details for a longer 
publication [TT| and now provide only the final results. 
So doing, in the limit t 00, s fixed, we obtain 



^ P>l/2 
1 P<l/2 



(5) 



which clearly points out the existence of a critical p — 
Pc = 1/2. The location of the critical point at p = 1/2 
is due to the fact that at this value, each particle on 



the boundary produces exactly one offspring. Below this 
point, the number of particles decreases exponentially 
with the number of past boundary particles, and above 
it it grows exponentially. The behavior of A(s) and a for 
the off-critical state, p ^ pc, is given by 



Ais) 



p> 1/2 



a 



1 + 13 p< 1/2 
(3 ]5>l/2 



(6) 



(7) 



The presence of the phase transition as we approach p = 
1 /2 from below is clearly manifest in the divergence of the 
coefficient A. The approach from above is not obvious 
from the large t behavior. What happens is that for 
p > 1/2, Pe first rises toward unity, as happens below 
the transition. However, at very large t, the behavior 
crosses over toward the power-law decay toward P^ . The 
details of this crossover will be presented in our longer 
presentation [TT]. 

In the critical state, p = Pc Eq. ([s]) remains valid and 
so P^ = 1 while the behavior of A{s) and a shows a 
transition as a function of /3. For a, we obtain 



1-/9 /3< 1/2 
13 f3>l/2 



(8) 



We can exhibit an analytic expression for A(s) at Pc only 
for 13 > 1/2: 

^ vr csc(/?.)r(l - b) 

^ ' r(i - 2/3)r(i + /3) ' ' ^ ' 

For (3 < 1/2, A{s) has to be calculated numerically in 
general. However, in the limit of large s, we have 



A{s) oc s 



W < 1/2). 



(10) 



The existence of a special /3 for the behavior at criticality 
is very non-trivial and it is especially interesting that the 
critical /3c is equal to 1/2, i.e the normal diffusion case. 
We can treat this critical f3c as a critical fractal dimension 
since 2(3 is just the fractal dimension of the diffusion. For 
the special case of /3 = 1/2 our ansatz, i.e. Eq. Q, does 
not work and one needs to treat this case specially [llj : 
the result is 



1-P£;(<,s)cx 



1 



log(s/i) 



s/t) 



1/2 



(/3 = l/2). (11) 



The logarithmic corrections in the behavior support our 
claims as to the critical nature of (3 ~ 1/2. For a two 
time quantity like the survival probability Ps{t,s) — 
1 — PE{t,s), one usually expects for a stationary pro- 
cess a dependence only on the time difference s; when 
the process is non-stationary this is generally not true. 
When the time dependence is that of a ratio of the two 
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times this is usually defined as aging, as in our case 
Ps{t,s) oc (s/ty-^. The aging behavior for /3 < 1/2 at 
criticality is a signature of the nonequlibrium phase tran- 
sition properties, which have been studied extensively in 
the context of DP and Contact Process [SJ [TS] , and 
the nonstationarity of the process usually obtained in 
glassy dynamics '1^. We must note that the obtained re- 
sult where the probability to observe at least one particle 
in time interval s is proportional to s^~^ (/3 < 1/2) is im- 
expected in light of the fact that for the off-critical state 
(p < Pc), this probability is proportional to s. 

Now we treat the limit s ^ t ^ 1, where we take first 
s — > oo. In this limit we can neglect the second term 
on the right hand side of Eq. (|3| and define PE{t,oo) = 
Psit) simply as the extinction probability of the process. 
We are now dealing with a single-time quantity and so 
no aging behavior could be obtained. For the solution 
we use the same ansatz as in Eq. Q, writing A{oo) ~ A. 
The equation for P|P remains the same as Eq. ([5| and 
one again obtains the same critical value for pj, = 1/2. 
In the off-critical state the results for a and A are 



and 



A{s) = 



p> 1/2 



1 

l-2p /? 



(12) 



(13) 



The transition in this limit is different then in the pre- 
viously discussed limit of t ^ s, as ^ diverges similarly 
near and thus describes the crossover time scale on 
which one would observe critical behavior. For the criti- 
cal state p = Pc we obtain 



PB(t)^l- 



/3/2 



(14) 



Thus for the survival probability of the process, Ps(t) = 
1 — Psit), we obtain the power- law behavior Ps{t) ~ 
t^^ with 6 = 13/2. The special properties of /3 = 1/2 
are not observed in the limit of s ^ t, and as already 
been mentioned no aging behavior could be observed for 
a one-time quantity. The noncomutativity of the limits 
t— >-cx)-s— ^cx)is similar in some sense to the non 
commutativity of the same limits present for correlation 
function behavior in glassy systems [1] and leads to such 
non-equilibrium property as ergodicity breaking [16l flT] 
Thus, in the case of normal diffusion, i.e. P = 1/2, 
for our modified model with unrestricted occupation of 
the origin, we have a phase transition at p = 1/2 with a 
critical exponent for the survival probability of (5 = 1/4, 
observed also in numerical simulations [2] . This is in 
contrast with the original model, with a larger critical p, 
but more importantly, a survival exponent of (5 ~ 0.17. 
The change in the phase transition point is clear, since in 
the original iV = 1 model, some children die upon their 
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FIG. 1: The survival probability vs. time for A'' = 1, 2, 3, 
and 4, together with power-law fits for the cases A'' = 1 and 
TV = 4. 



attempted return to the origin and so the critical p has 
to be larger than 1/2 to compensate. The change in ex- 
ponent is not unexpected. For directed percolation, the 
critical exponents in a model of unrestricted occupancy 
are the mean-field exponents, which differ from the ex- 
ponents in the restricted occupancy version[18 . 

In the case of directed percolation, for any finite N 
the critical exponents are those of the = 1 model. For 
large A^, there is a crossover in the scaling between the 
mean-field scaling of the N = oo model and the directed 
percolation exponents [18J . Based on this analogy, one 
would expect a similar behavior in the present case, with 
the critical behavior for any finite A^ being identical to 
the A^ = 1 model. 

To test this hypothesis, we have measured the survival 
exponent S at criticality for various A^. The results are 
shown in Fig. [T] We see that, contrary to our naive 
expectation, there is a different exponent for each N. The 
measured best fit exponents are, for example, S = 0.160 
for A^ = 1, (5 = 0.171 ioT N = 2,6 = 0.180 for A^ = 3 and 
6 = 0.192 for A^ = 4. There is no sign of any crossover 
behavior. One test for this is shown in Fig. [2j where 
we show Ps{t)t^ , in a linear(?/)-log(a;) scale. We see that 
the graphs all exhibit the same characteristic oscillation 
in ln(x) [19] that makes accurate determination of an 
exponent so difficult, but they show no secular trend. 

Varying A^ is not the only way to interpolate between 
the two models and thereby vary 5. We can instead intro- 
duce a parameter A which controls the lifetime of a par- 
ticle on the active site after reproduction (in the original 
model, A = 1). As A decreases, the interference between 
different particles is reduced, and the theory approaches 
the noninterference, Galton- Watson, limit. We can also 
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FIG. 2: The survival probability, , multiplied by e"*"* as 
a function of time for A'' = 1, 2 and 3. The values of 5jv are 
as noted in the text. 

increase A beyond one, and make the interference ef- 
fect stronger, thereby reducing the exponent below the 
original BH value. Actually the precise dimensionless pa- 
rameter which controls the exponent S is cx -DA, 
which presents a surprising dependence of the critical ex- 
ponent 6 on the diffusion coefficient (D) in the bulk. This 
is demonstrated in Fig. [] where the survival exponent is 
plotted as a function of A. We see that as opposed to 
the discrete parameter N, now the exponent varies con- 
tinuously, decreasing with DA. 

Finally we note that the relevance of the limit for 
which we managed to solve the model analytically is 
much broader when one considers higher dimension sys- 
tems. We have studied herein a one dimensional system 
where the boundary is zero dimensional, while if one con- 
siders a d dimensional system where the boundary is of 
dimension d — 1, a line or a membrane, in such a case 
even if DA does not approach zero, the probability that 
two particles would try to occupy the same site (since 
there an infinite number of sites on boundary) during 
the time DA is small and this would lead to a larger 
applicability of our analytical solution. Another type of 
system where our results could be applied is systems with 
anomalous diffusion behavior described by the Continu- 
ous Time Random- Walk model p3] , for such systems the 
results stays valid even if the size of the system is finite 
as for a finite system with regular diffusion our results 
are modified for times long enough [TT] . 

In summary, we have extended to finite active site oc- 
cupancy N and solved analytically the N ^ oo limit of 
the Barato-Hinrichsen model, yielding the scaling expo- 
nents of the extinction transition in this limit. We have 
also show that below a critical fractal dimension for the 



FIG. 3: The measured survival exponent (5 as a function of 
A, the lifetime of the particle on the active site after it gives 
birth. The simulations were performed for measurement time 
t ~ 10^^ 



diffusion process in the bulk, there is aging behavior at 
the transition. In addition, we have seen that the expo- 
nents are in fact not universal, varying with N or the 
scaled lifetime of the particles on the active site. This is 
of course very different from the situation that obtains in 
the superficially similar contact process, which exhibits 
universal Directed Percolation scaling. 
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